The influence of the hydrodynamic drag from an accretion torus on extreme mass-ratio inspirals 
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We have studied extreme mass-ratio inspirals (EMRIs) in spacetimes containing a rotating black hole and 
a non self-gravitating torus with a constant distribution of specific angular momentum. We have found that 
the dissipative effect of the hydrodynamic drag exerted by the torus on the satellite is much smaller than the 
corresponding one due to radiation reaction, for systems such as those generically expected in Active Galactic 
Nuclei and at distances from the central supermassive black hole (SMBH) which can be probed with LISA. 
However, given the uncertainty on the parameters of these systems, namely on the masses of the SMBH and of 
the torus, as well as on its size, there exist configurations in which the effect of the hydrodynamic drag on the 
orbital evolution can be comparable to the radiation-reaction one in phases of the inspiral which are detectable 
by LISA. This is the case, for instance, for a 10''A/q SMBH surrounded by a corotating torus of comparable 
mass and with radius of 10^ — 10'' gravitational radii, or for a lO^i\f0 SMBH surrounded by a corotating 
W^Mq torus with radius of 10^ gravitational radii. Should these conditions be met in astrophysical systems, 
EMRI-gravitational waves could provide a characteristic signature of the presence of the torus. In fact, while 
radiation reaction always increases the inclination of the orbit with respect to the equatorial plane (i.e., orbits 
evolve towards the equatorial retrograde configuration), the hydrodynamic drag from a torus corotating with the 
SMBH always decreases it (i.e., orbits evolve towards the equatorial prograde configuration). However, even 
when initially dominating over radiation reaction, the influence of the hydrodynamic drag decays very rapidly 
as the satellite moves into the very strong-field region of the SMBH (i.e., p < 5A/), thus allowing one to use 
pure-Kerr templates for the last part of the inspiral. Although our results have been obtained for a specific class 
of tori, we argue that they will be qualitatively valid also for more generic distributions of the specific angular 
momentum. 
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I. INTRODUCTION 

One of the most exciting prospects opened up by the sched- 
uled launch of the space-based gravitational-wave detector 
LISA lU will be the possibility of mapping accurately the 
spacetime of the supermassive black holes (SMBHs) which 
are believed to reside in the center of galaxies [2]. Among 
the best candidate sources for this detector there are Extreme 
Mass Ratio Inspirals (EMRIs), i.e. stellar-mass black holes 
(771 « 1 — IQAIq) or compact objects orbiting around the 
SMBH and slowly inspiraling due to the loss of energy and an- 
gular momentum via gravitational waves (radiation reaction). 
In order for the signal to fall within the sensitivity band of 
LISA, the SMBH must have a mass M « 10^ - IO^A/q, i.e., 
the low end of the SMBH mass function. 

It is currently expected that a number of such events ranging 
from tens to perhaps one thousand could be measured every 
year yjl, but since they will have small signal-to-noise ratios, 
their detection and subsequent parameter extraction will re- 
quire the use of matched-filtering techniques. These basically 
consist of cross-correlating the incoming gravitational-wave 
signal with a bank of theoretical templates representing the 
expected signal as a function of the parameters of the source. 
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This will not only allow one to detect the source, but also 
to extract its properties. For instance, the accurate mod- 
elling of the motion of a satellite in a Kerr spacetime will 
allow one to measure the spin and the mass of the SMBH. 
Although producing these pure-Kerr templates has proved to 
be a formidable task, particularly because of the difficulty of 
treating rigorously the effect of radiation reaction (see Ref. fl 
for a detailed review), considerable effort has gone into trying 
to include the effects of a deviation from the Kerr geometry. 
These attempts are motivated by the fact that possible "ex- 
otic" alternatives to SMBHs have been proposed (e.g., boson 
stars p"], fermion balls [6] and gravastars [7]), although the 
presence of these objects would require to modify radically 
the mechanism with which galaxies are expected to form. On 
the other hand, non-pure Kerr templates might allow one to 
really map the spacetimes of SMBHs and to test, experimen- 
tally, the Kerr solution. 

Different approaches to this problem have been considered 
in the literature. EMRIs in a spacetime having arbitrary grav- 
itational multipoles should be considered in order to maintain 
full generality |8], but this method does not work very well in 
practice and would only apply to vacuum spacetimes. For this 
reason, alternative approaches have been proposed and range 
from EMRIs around non-rotating boson stars [9], to EMRIs in 
bumpy black-hole spacetimes 1 10] (i.e., spacetimes which are 
almost Schwarzschild and require naked singularities or ex- 
otic matter) or in quasi-Kerr spacetimes li 1 lil (i.e., spacetimes. 



consisting of Kerr plus a small quadrupole moment). 

Interestingly, none of these methods is suitable for taking 
into account the effect of the matter which is certainly present 
in galactic centers. SMBHs can indeed be surrounded by stel- 
lar disks (as in the case of the Galactic center I.1Z1 ) or, as in the 
case of Active Galactic Nuclei (AGNs) 11311 . in which we are 
most interested, by accretion disks of gas and dust which can 
be even as massive as the SMBH 1141]. While the gravitational 
attraction of a disk can have important effects on EMRIs if 
this disk is very massive and close to the SMBH |15], an as- 
trophysically realistic accretion disk can influence an EMRI 
only if the satellite crosses it, thus experiencing a "hydrody- 
namic" drag force. 

This drag consists of two parts. The first one is due to 
the accretion of matter onto the satellite black hole (this was 
studied analytically by Bondi & Hoyle [,16,1 an d subsequently 
confirmed through numerical calculations Ill7l[l8l[l9ll ). This 
transfers energy and momentum from the disk to the satellite, 
giving rise to a short-range interaction. The second one is 
instead due the gravitational deflection of the material which 
is not accreted, which is therefore far from the satellite, but 
which can nevertheless transfer momentum to it. This long- 
range interaction can also be thought of as arising from the 
gravitational pull of the satellite by its own gravitationally- 
induced wake (i.e., the density perturbations that the satel- 
lite excites, by gravitational interaction, in the medium), and 
is often referred to as "dynamical friction". This effect was 
first studied in a collisionless medium by Chandrasekar ll20ll . 
but acts also for a satellite moving in a collisional fluid 

([nllllllllllllilll. 

The effect of this disk- satellite interaction on EMRIs has 
been studied by different authors for a number of disk mod- 
els. More specifically, in a series of papers Karas, Subr and 
VokrouhUcky considered the interaction between stellar satel- 
lites and thin disks |J2^ HHll. In Ref. Ill, in particular, 
Subr and Karas found that the effect of the star-disk inter- 
action on EMRIs dominates over radiation reaction for thin 
disks, both for non-equatorial orbits crossing the disk only 
twice per revolution and for equatorial orbits embedded in the 
disk. The only exceptions to these conclusions come if the 
satellite is very compact (a neutron star or a black hole) or the 
disk has a low density {e.g., in the region close to the SMBH if 
the flow becomes advection-dominated). These results agree 
with those found by Narayan 1291 , who focused on Advection 
Dominated Accretions Flows (ADAFs), which were believed 
to describe accretion onto "normal" galactic nuclei (i.e., ones 
much dimmer than AGNs)'. Overall, he found that for com- 
pact objects and white dwarfs the effect of the hydrodynamic 
drag is negligible with respect to radiation reaction, whereas 
it is not negligible for main sequence and giant stars. More 
recently. Levin |31] has proposed a scenario in which mas- 



' Accretion onto "normal" galactic nuclei is now believed to be better de- 
scribed by Advection Dominated Inflow Outflow Solutions (ADIOS) 13011 . 
However, this is not expected to change significantly Narayan's results 
since ADIOS's, like ADAFs, have very low densities in the vicinity of the 
SMBH. 



sive stars form in a thin accretion disk in an AGN, ultimately 
producing stellar-mass black holes embedded in the disk. The 
small black holes are then dragged towards the (non-rotating) 
SMBH, but if this is accreting at a rate comparable to the 
Eddington limit, the drag from the accreting gas will not af- 
fect the final part of the inspiral (i.e., at radii smaller than 10 
Schwarzschild radii) significantly. Finally, Chakrabarti ll32ll . 
studied instead the orbital evolution of a satellite black hole 
on a circular equatorial orbit embedded in a disk with a non- 
Keplerian distribution of angular momentum, and found that 
the exchange of angular momentum between the disk and the 
satellite can lead to significant orbital modifications. 

All of these studies have been been carried out within a 
Newtonian or pseudo-Newtonian description of gravity (with 
the partial exception of Ref. |26], in which the orbits are Kerr 
geodesies, but the disk model and the hydrodynamic drag is 
not relativistic). In this paper, instead, we provide a first rel- 
ativistic treatment for satellite black holes moving on generic 
orbits around a rotating SMBH surrounded by a thick disk 
(i.e., a torus). We consider the torus to have constant specific 
angular momentum and neglect its self-gravity (i.e., we con- 
sider the metric to be pure Kerr). Under these assumptions, an 
analytical solution exists for this system ll33U34ll . This con- 
figuration can be proved to be marginally stable with respect 
to axisymmetric perturbations |35] (i.e., if perturbed, such a 
torus can accrete onto the SMBH), and is expected to be a 
good approximation at least for the inner parts of the accre- 
tion flow iM HI ■ 

We have found that for a system composed of a SMBH with 
mass AI — 10^ Mq and a torus with mass A/t < AI and outer 
radius rout = lO^M, the effect of the hydrodynamic drag 
on the motion of the satellite black hole is much smaller than 
radiation reaction at those distances from the SMBH which 
can be probed with LISA (i.e., - lOAf for M = lO^Af©). 
Although these values for A/, A/t and rout are plausible for 
AGNs, an overall uncertainty is still present and has moti- 
vated an investigation also for different masses and sizes of 
the torus. In this way we have found that the effect of the 
torus can be important in the early part of the inspiral and 
that it could leave an observable imprint in the gravitational 
waveforms detected by LISA, if the radius of the torus is de- 
creased to rout = 10"^ — 10'' Af or, even for rout = lO^Af and 
A/t < M, if M = 10^ Mrs,. I" this latter case, in fact, LISA 
could detect an EMRI event at distances as large as r ^^ 45A/ 
from the SMBH, although the event needs to be sufficiently 
close to us because the amplitude of the gravitational-wave 
signal decreases as A//r. 

In addition, if non-negligible, the effect of the hydrody- 
namic drag would have a distinctive signature on the wave- 
forms. Radiation reaction, in fact, always increases the in- 
clination of the orbit with respect to the equatorial plane 
(i.e., orbits evolve towards the equatorial retrograde config- 
uration) IBql . The hydrodynamic drag from a torus corotat- 
ing with the SMBH, on the other hand, always decreases this 
angle (i.e., orbits evolve towards the equatorial prograde con- 
figuration). Should such a behavior be observed in the data, 
it would provide a strong qualitative signature of the presence 
of the torus. However, it is important to point out that even for 



those configurations in which the hydrodynamic drag plays a 
major role, this is restricted to the initial part of the inspiral 
detectable by LISA, whereas its effect rapidly vanishes in the 
very strong-field region of the SMBH (i.e., p < 5M). As a 
result, the pure-Kerr templates would provide a faithful de- 
scription of the last part of the inspiral even in these cases. 

The rest of the paper is organized as follows. In Sec. |II] 
we review the equilibrium solutions that we used for the or- 
biting torus. In Sec. IIII Al we present the equations governing 
the interaction between the satellite black hole and the torus, 
while in Sect. IIII B] we apply the adiabatic approximation to 
the hydrodynamic drag. Results are then discussed in Sec. 
IIV Al for equatorial circular orbits and in Sec. lIVBI for generic 
(inclined and eccentric) orbits. Finally, the conclusions are 
drawn in Sec.|V] Throughout this paper we use units in which 
G = c=l. 



II. MODELLING THE TORUS 

The properties of non self-gravitating, stationary, axisym- 
metric and plane-symmetric toroidal fluid configurations in 
Kerr spacetimes are well-known in astrophysics but are less 
well known within the community working on EMRIs. Be- 
cause of this, in this Section we briefly review the basic facts, 
referring the interested reader to Refs. iH [M Hzl [M HHl f or 
additional information. 

Let us consider a perfect fluid with 4- velocity it^""^, which 
is described by the stress-energy tensor 



T'^'' = (P + P)"(^uid<uid+P5'^'' 



(1) 



where p, pq, p and h= {p + p)/po are the pressure, rest-mass 
density, energy density and specific enthalpy of the fluid. In 
what follows we will model the fluid with a polytropic equa- 
tion of state p ~ kpq — poe{r — 1), where e ~ p/po ^ 1 
is the internal energy per unit rest-mass, and k and F are the 
polytropic constant and index, respectively. Because we are 
neglecting the self-gravity of the fluid, we can also consider 
g as given by the Kerr metric in Boyer-Lindquist coordinates, 
which reads li4Ql 



ds^ = - 1 - 
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where the second equals sign underlines that the vector and 
the 1-form are each the dual of the other. Here, Vt = 



"fluid /"fluid is the angular velocity, A 



.fluid 



is called the 



is the energy per unit mass as 



^fluid 



— „, fluid 



^yfluid jg j.jjg specific an- 



redshift factor, U = 
measured at infinity and 
gular momentum as measured at infinity {i.e., the angular mo- 
mentum per unit energy as measured at infinity). Note that t is 
conserved for stationary axisymmetric flows, as can be easily 
shown using Euler's equation. The specific angular momen- 
tum and the angular velocity are trivially related by 






gtt + gt<t,^ 



(5) 



while the normalization condition u^^^'^ ■ u 



fluid 



U 



A = 



I gtd^ + '2gt4,e + 500 
/ -1 

9tt + '2gt4>^ + g4,ci>^^ 
1 



AU 



1 -Qi 
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(6) 

(7) 
(8) 



where w'^ ~ g^i — gu g^^ = A sin^ 9. Note that in this paper 
we will always consider i > {) (torus rotating in the positive 
(/)-direction), while we will allow the spin parameter a of the 
black hole to be either positive (black hole corotating with the 
torus) or negative (black hole counter-rotating with respect to 
the torus). 

To calculate the structure of the torus, we need to use Eu- 
ler's equation, which in its general form reads 



(.<" 



^fluid "fluid 



)duP 



"-fluid 



P 



(9) 



where afj^^j^^ is the 4-acceleration of the fluid. In particular, if 
the pressure is assumed to depend only on r and 9 and if the 
equation of state is barotropic [i.e., \f p — p(p)]^, from Eq. 
(|9]l one easily gets that the 4-acceleration can be expressed as 
the gradient of a scalar potential W{p): 



fluid 



= d^w , w{p) = - 



dp' 



p' + p{p') 



(10) 



On the other hand, from the definition of 4-acceleration 
«uid = "fluidV,.M(^uid)' Eqs. ®, © and ®, and the Kifling 
equation V(^^^) = for ^ = d/dt and ^ = 9/90, one easily 
gets 



where 



a^ cos^ ( 



A; 



2Mr 



(3) 



The fluid is assumed to be in circular non-geodesic motion 
with 4-velocity 



.fluid 
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dt 
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d 
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In particular, taking the derivative of this equation, anti- 
symmetrizing and using the trivial fact that i9[^^]VF = 



^ This is of course the case for a polytropic equation of state, because p = 

p/(r - 1) + (p//.)i/r. 



d[nu]^ — d[^iy]U = 0, we obtain that d[^^dy]li — 0. This 
impHes Vfi c>c Vt and thus that I and Vl have the same con- 
tour levels {i.e., Q. = f2(£)]. Using this fact, we can then write 
Eq. (fTTT l in an integral form: 



W -Wo 



= In [/ - hi U„ 



dp' 



P' 



tau 



■ Pip') 

n{e')de 
^ 1 - n{£')£' 



(12) 



where Wout and £out are the potential and specific angular 
momentum at the outer edge of the torus. ^ 

In the case of a torus with constant specific angular momen- 
tum [i.e., i{r, 9) — constant], Eq. (fTSl i provides an analytical 
solution, because once i has been fixed the integral on the 
right-hand side is zero and Eq. (|6]l gives an analytical expres- 
sion for U: 



W-Wo. 



dp' 



p' + p{p') 



\nU -\nUo 



(13) 



Note that if one requires that W ^ Q when r -^ +cx) {i.e., 
W = {) for an equipotential surface closing at infinity), this 
equation gives W = \iiU: W > i) then corresponds to open 
equipotential surfaces, while W < Q corresponds to closed 
equipotential surfaces. Interestingly, the potential well can 
present a minimum and a saddle point. Because of the plane- 
symmetry, these points are located in the equatorial plane, 
thus corresponding to local extremes of W{r, — 7r/2), and 
mark two important positions: respectively, the center of the 
torus {i.e., the point where the density reaches its maximum) 
and its cusp {i.e., the mass-shedding point). Noticeably, these 
points are located at the radii where the specific angular mo- 
mentum of the torus, i, coincides with that of the geodesic cir- 
cular equatorial orbit (the "Keplerian" orbit) corotating with 
the torus. 



iK{r,a) 



{r-2M)y/r/M + a 



(14) 



This immediately follows from the fact that at the extremes 
of the function W one has df^W = 0, which leads, through 
Eq. ([Tol l, to a^^jjj — (in other words, at the center and at 
the cusp the pressure gradients are zero and only gravitational 
forces act). 

In this paper we will indeed consider constant-^ tori. A de- 
tailed classification of these models depen ding upo n the val- 
ues of £ and of Wout can be found in Refs. |l33ll34ll37ll . Here 
we simply recall that in order to have a closed equipotential 
surface with a cusp, one needs to have a value of £ between 
the specific angular momenta ^,„s and f ,nb of the marginally 
stable and marginally bound equatorial geodesic {i.e., "Ke- 
plerian") orbits corotating with the torus. This can be eas- 
ily understood by noting, from Eq. ( fTSl ), that the potential 



^ Of course, Wout and iout can be replaced by the values of W and I at the 
inner edge of the torus if this is present. 



W{r, — 7r/2) is simply the effective potential describing 
the equatorial motion of a test particle around a Kerr black 
hole. As such, £ms and £mb can be calculated easily using Eq. 
(fT4l) and the formulas for the radii of the marginally stable and 
marginally bound circular equatorial orbits in Kerr rotating in 
the positive (/)-direction {i.e., corotating with the torus): 

^ms = iK{rms) , ^mb = ^x(?'mb) , (15) 

r^jM - 3 + ^2 - sign(a)V(3-^i)(3 + Zi + 2Z2), 



(16) 
(17) 



r^b/M = 2 - a + 2^1 - a , 

Zi = 1 + (1 - ^2)1/3 [(1 + 5)1/3 + (1 _ 5)1/3] ^ (18) 
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(19) 



where a ~ a/M. 

In order to pick up a particular solution having both an inner 
and an outer radius, one needs also to choose a negative value 
for the "potential barrier" at the inner edge of the torus. 



AW = Ww 



Wc, 



Wo, 



W^cusp < . 



(20) 



If AW < 0, the inner radius of the torus is larger than the 
radius at which the cusp occurs (rjn > rcusp), while if the 
potential barrier AW reduces to zero, the torus exactly fills 
its outermost closed equipotential surface and ri„ = 7-cusp < 
Tins- Note that because of the considerations that we have 
made above about the value of £, for constant-^ tori we have 
^-cusp > r-nib (with Tcusp = r^b Only if £ = £mb) and 
reenter > ^-jiis (with Tccntor = ^-ms Only if £ = £ms)- If instead 
AW > 0, the fluid overflows the outermost closed equipo- 
tential surface and mass transfer is possible at the cusp: for a 
polytropic equation of state, the accretion rate can be shown 
to be M ex AM^i'/^r-i). 

The integral Euler equation for constant-^ tori [Eq. (fTsl ll 
further simplifies if the equation of state is polytropic, because 
in this case 



dp' 



P' + P{P') 



In- 



ho 



(21) 



where /lout is the specific enthalpy at the outer edge of the 
torus. Since for a polytropic equation of state the enthalpy is 
given by 



1 



-kpI ^ 



(22) 



it is clear that /lout — 1 (because p = po = at the outer edge 
of the torus), and Eqs. ( fTlT l and ( 1211 1 give 



Po(r-,0) = 



r-1 [e 



Waut-W(r,e) 



1] 



i/(r-i) 



(23) 



Once the rest-mass distribution is known, the total rest mass 
of the torus is given by 



Mt,Q = 



/—gu d X , 



(24) 



where ^/— g = EsinS and d^x = dr dO d(j) is the coordinate 
3-volume element, while the mass-energy reads 



27r 




po>0 \94>4> + '^9t4,^ + 9tti'^ 
{r'^ +0^ cos^e) sine drde . 



(25) 



Clearly, the smaller the ratio between the mass of the torus and 
that of the SMBH, the better the approximation of neglecting 
the self-gravity of the torus. 



III. MODELLING THE ORBITAL MOTION 

This Section is dedicated to the discussion of the hydro- 
dynamic drag on the satellite black hole. Although the two 
aspects are closely inter-related, we first discuss the equations 
governing the interaction between the satellite black hole and 
the torus and then describe their use in the calculation of the 
changes of the orbital parameters within the adiabatic approx- 
imation. 



A. The hydrodynamic drag 

As already mentioned in Sec. |I] the hydrodynamic drag act- 
ing on the satellite black hole can be written as the sum of a 
short-range part, due to accretion, and a long-range part, due 
to the deflection of the matter which is not accreted or, equiv- 
alently, to the gravitational interaction of the satellite with the 
density perturbations gravitationally induced by its own pres- 



dr 



dp^' 
"d7 



dpi" 
"d7 



dcfi 



(26) 



where r is the proper time of the satellite. 

Accretion onto a moving black hole was studied analyti- 
cally in a Newtonian framework by Bondi & Hoyle 1161 . who 
found the rest-mass accretion rate to be 



dmo 
~d^ 



(w2 +tl2)3/2 



(27) 



where m is the mass of the black hole, v and Vg are respec- 
tively the velocity of the black hole with respect to the fluid 
and the sound velocity, and A is a dimensionless constant of 
the order of unity, which for a fluid with polytropic equation 
of state and polytropic index F has the value |41] 



j^x (r+i)/[2(r-i)] 
2J 



5-3r 



-(5-3r)/[2(r-i)] 



(28) 



Subsequent numerical work iITtI [Tsl \l% treated instead the 
problem of accretion in full General Relativity, and showed 
that Eq. ( |27] i. with A given by Eq. ( |28] l. is correct provided 
that it is multiplied by a factor ^^ 5 - 25 when v and Vg become 



relativistic (cf. Table 3 of Ref. lUSll ). However, because a 
fit for this correction factor is, to the best of our knowledge, 
not yet available, and the published data is not sufficient for 
producing one, we use the Bondi accretion rate [Eqs. ( l27b and 
(|28]|1. bearing in mind that it could slightly underestimate the 
drag at relativistic velocities v and Vg^ Once the accretion 
rate is known, the short-range part of the drag reads iflTIl 



dp^' 
"d7 






(29) 



where we recall that h is the specific enthalpy of the fluid. 
Note that this equation basically follows from the conserva- 
tion of the total 4-momentum of the satellite and the fluid. 

The long-range drag is instead more complicated. The 
gravitational interaction of a body with the density perturba- 
tions that it excites gravitationally in the surrounding medium 
was first studied by Chandrasekhar |20] in the case of a colli- 
sionless fluid, and is also known as "dynamical friction". Al- 
though less well recognized, dynamical friction acts also for a 
body moving in a coUisional medium lll7ll2l1. 22, 23, 24, 25]. 
In particular, a satellite moving on a circular planar orbit {e.g., 
a circular orbit around a Schwarzschild black hole or a circu- 
lar equatorial orbit around a Kerr black hole) experiences a 
drag in the tangential direction Ill7ll2lll22ll23ll and one in the 
radial direction [24|l: 
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defl 
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dcfl 
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dofl 
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(30) 



where cr is a unit spacelike vector orthogonal to Msat and 
pointing in the direction of the motion of the fluid, 



"fluid - T^sat 



Vt^ 



1 



(31) 



(the Lorentz factor 7 = —"fluid • "sat encodes the relative 
motion of the satellite with respect to the fluid of the torus), 
and 



M^'^'Usat -Or<T 



d/dr 



[5„-(u-t)V(72-l)]i/2' 



(32) 



is a unit spacelike vector, orthogonal to both Wgat and a and 
poiting in the radial direction. In particular, the tangential and 
radial drags are given by 1124,12511 
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dp 
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dcfl 
rad 



dofl 



tang 1 



Ait{p + p)m^j^{l 



.,2\2 



■ -'rad ; 



(33) 



(34) 



'* As we will see in section IIVI v and Vs can become relativistic only for 
orbits counter-rotating with respect the torus and very close to the SMBH. 
For these orbits the dominant part of the hydrodynamic drag is the long- 
range one, and the relativistic correction factor to the Bondi accretion rate 
(which is roughly 5 — 10 for these orbits, as can be seen comparing the 
middle panel of Fig. [T]with Table 3 of Ref. llSll ^ does not change this 
conclusion. 



where /tang and /lad are complicated integrals. Fits to the 
numerically-com pute d steady-state'' values for these integrals 
are given in Ref . ll24ll : 



0.7706 In 



ir 



1+M 



0004-0. 9185M 



1.4703A^, 



forAl < 1.0, 



'tang 



and 



ln[330(r/r„,i„)(A^ - 0.71)^-72_yy^-9-58] 
for 1.0 <M <4:A, 

M(r/r„,i„)/(0.11M + 1.65)], 
^ for M>4A, 



' A^2 ^0 3-51-^-4.22^ forX<l.l, 

- ^ I 0.5 In [9.337W2(7W2 - 0.95)], 
'■^^ ^ ] for 1.1 <7W <4.4, 

.0.3 7^2^ forX>4.4, 



(35) 



(36) 



where r is the radius of the circular orbit, rmin ~ 2m(l + 
t|2)/t|2 is the capture impact parameter of the satellite black- 
hole, while M = v/vs is the Mach number. 

These fits are valid for r ^ rmin and are accurate within 
4% for M < 4.4 and within 16% for M > 4.4. However, the 
fit for /tang does not go to zero when M goes to zero, while 
^lad goes to zero only as Ai"^ in this limit: these behaviors 
would give a non-zero radial drag and a diverging tangential 
drag for w — > [c/ Eqs. ( |33] ) and (l34bl. This is clearly a spuri- 
ous behavior: dynamical friction must vanish for v — 0, since 
in this case the pattern of the density perturbations is spheri- 
cally symmetric around the body (as there is no preferred di- 
rection). However, as we will see in Sec. HUB I the effect of 
the radial drag vanishes if one uses the adiabatic approxima- 
tion (as it is usually done in EMRI-studies ll42ll43ll44ll45ll ). 
and therefore this artifact of the fit ( |36] | cannot cause any harm 
in our numerical code. This is instead not the case for the 
tangential drag: in order to eliminate its spurious divergence, 
we have approximated /tang with its straight-line functional 
form at low Mach numbers. Since the dynamical friction drag 
for straight-line subsonic motion is given by Eq. ( [33] ) with 

-'tang 



1/2 ln[(l + M)I{1 -M)]-Mk X 73 + 7W 75, 



we can assume that /tang is given, for M. < 0.1, by 



^tang 



0.9563 



M^ 



7W5 
~ 



(37) 



Note that although Eq. ( |30] | is strictly valid only for circular 
planar motion {i.e., in the case of a Kerr spacetime, for circu- 
lar equatorial orbits), we expect it to be a good approxima- 
tion also for generic orbits around a Kerr black hole. Indeed, 
thanks to the choice of the unit vectors <j and x, Eq. (|30] | 
gives a tangential drag parallel the direction of the flow and 
a drag in the radial direction perpendicular to the direction of 
the flow. Both of these components are expected to be present 
also for generic orbits. In particular, the tangential drag should 
be given approximately by Eqs. ( [33] ) and ( [37] l if the radius r 
appearing in Eq. ( [35] ) is replaced by the semi-latus rectum 
p of the orbit [see Eq. ([52] | for the definition of this quan- 
tity].^ Although this prescription is not exact, the results of 
Ref. II24I1 suggest that the relevant lengthscale in the Coulomb 
logarithm appearing in the second and third lines of eq. ( [35] ) 
should be one characterizing the orbit, rather than the size of 
the medium, as commonly assumed in most of the works on 
dynamical friction predating Refs. [22,^?] (see the introduc- 
tion of Ref. 1 24] and references therein for more details about 
this point). Of course, this lengthscale could be different from 
the semi-latus rectum of the orbit, but different choices for it 
would have only a slight impact on the results because of the 
logarithmic dependence. 

The extrapolation of the radial drag given by Eqs. ( [34] | and 
( [36b from circular planar to generic orbits is instead a bit more 
problematic, although one expects it to be a good approxima- 
tion at least for orbits with small eccentricities and small in- 
clinations with respect to the equatorial plane. At any rate, as 
we have mentioned earlier, in Sec. IIIIBI we will show that the 
effect of this radial drag on the orbital evolution averages to 
zero when adopting the adiabatic approximation. (Note that 
this agrees with Ref. |24], which found that the effect of the 
radial drag on the orbital evolution was subdominant with re- 
spect to that of the tangential drag.) Nevertheless, a non-zero 
effect may still be present in cases in which the adiabatic ap- 
proximation is not valid (i.e. if the hydrodynamic drag acts 
on a timescale comparable to the orbital period), or possibly 
even in the adiabatic approximation if more rigorous expres- 
sions for the radial drag should be derived in the future. 

The rate of change of the mass of the satellite with respect 
to the coordinate time t follows immediately from p^^^. p^J^^ — 
—m?: denoting the derivative with respect to t with an over- 
dot, we have 



^.sat 7 u 

uLt dr 



wLf dr 



hj drriQ 
uLf dr 



(38) 



It is well-known Il46ll that Kerr geodesies can be labeled, up 



where the factor 0.9563 is introduced to match the above fit at 
M = 0.1. 



^ Fortunately, the steady-state values for these integrals are reached over 
timescales which are comparable with either the sound crossing-time r/vs , 
r being the radius of the circular orbit, or with the orbital period. 



* Note that the tangential drag given by Eqs. (33), (35) and (37) is approx- 
imately connect also for straight-line motion, if r replaced in Eq. (35) by 
vt - t being the time for which the satellite has been active |22, 25] - as 
long as vt is smaller than the size of the medium, and by a cutoff-length 
of order of the size of the medium at later times. To see this, compare 
Eqs. (35) and (37) to the functional form of /tang for straight-line mo- 
tion, which is /tang = 1/2 ln[(l + A4)/(l — A^)] — M for subsonic 
motion and /tang = l/21n(l — 1/A4^) + ln(i)t/rniin) for supersonic 
motion El [111. 



to initial conditions, by three constants of motion, the dimen- 
sionless energy E and the angular momentum L^ as measured 
by an observer at infinity. 



L,=<7Af, 



E ^ 



and the dimensionless Carter constant Il46ll Q, 



(39) 



Q= M-j + a'^ cos^ 9(1 - E'^) + cot^ ell (40) 

where a — a/M . We will now derive expression for the rates 
of change of these quantities. 

To this purpose, let us first introduce the tetrad (Mgat, ^i = 
(T, 62 = X, 63) based in the position of the satellite and write 
the change in the 4-velocity due to accretion and deflection of 
the flow as 



-^uL 



(t) M 
-'sat'^^sat 



^sat^(i) 



(41) 



where Su^^^^ and Sujl^^ are the components with respect to the 



,W„,(J) 



teti-ad. In particular, perturbing — (Ugat)^ + '^y ^sat^sat = ^1 
to first order one easily gets -u^ll^u^ll + ^y^sat'^u^^i = 0, 
and using then the fact that Ug!j^ = to zeroth order, one 

obtains ^w^^t = 0- Using now Ju^^t = '^p^at/"^. e^i) ■ u^at = 
and e(j) • cq) — Sij (i — 1, 2, 3), Eq. (|4TI ) becomes 



The evolution of Q therefore foUows from Eq. (|40] | 



sat ' 



mo h 



m + ^L 

dE dL, 



^tang 



du 



^i + ^z. 

dE dL, 



771 -^/t^ 

dQ 



sat 
Ug 



dQ Pll^ 



■7 



TTlQ h 



dul 



• tang 
PdcR 



Xe 



my^ 



dul 



Xe, 



(47) 



where the partial derivatives are meant to be calculated with 
Eq. (|40l ). Note that the first and the second term of the first 
line cancel out because Q is conserved for geodesic motion. 

A useful alternative form for the evolution rate of Q can be 
obtained by rewriting Eq. ( |40l ) using the normalization condi- 



Q = A-^ E{f^ +d.^)-aL 



{L - aEf 



A«*) , (48) 



where f = r/M and A = A/Al"^. Proceeding as above and 
in particular using the fact that 



-^sat ~ 



Smo h 



SniQ h 



"%uid '^'^ 



X„tang 
"^dcfl 



+ 



^„tang 
"Pdofl 



my^ 



("fluid 



7<at) 



^,u 



^ u 



(42) 



Using now Eqs. dUl, (|39] | and (|42] |. we immediately obtain 



E 

E 




I. 



niQ h 



■ tang 
Pdcfi 



(. 



eu 



i/T^T / \ml. 



■rad 



•rad 
/'dofl 



(43) 



mML 



^X'P 



(44) 



In order to calculate instead the rate of change of the dimen- 
sionless Carter constant Q, let us note that from Eq. (|42] | it 
follows that the variation of ug in a short time interval St due 
to accretion and deflection of the flow is 



Sur = 



rfiQ h 



-tang 
Pdofi 



777^7^ 



1 



•rad 

uT' + ^^Xs 

777 



St. 



We can then write iig^^ as the sum of a term coming from the 
gravitational evolution (i.e., the geodesic equation) and one 
coming from collisions with the surrounding gas: 



M ,,sat^i/ 
sat 



i iiy'^H ^ 



777o h 



m 



^tang 
T^dcfl 



7)7, a/t^ 



■rad 

m 



xe- (46) 



^ ri/"/j ■^' 



sat^i/ 
sat 



7 



TTlQ h 



Pdcfl 



^ m 777^72 - ly 
[from Eqs. (|4|i and (l42]il. one easily gets 
dQ 



^; — r^ ^sat 

or 
9Q / 



777o /l 



dE 



•rad 
Pdcfl . 

7)7, 






(49) 



dui^* 



+ 



,^tang 
T^dcfl 



9Q S 



dE dL^ 



777 YT^ 



TTlQ h 



dul^* 



.tang 
Pdcfl 



Xr 



myj^ — 1 



5Q pS 



Xr 



(50) 



cJzi^'** 7)7 

where the partial derivatives are now calculated with Eq. ( |48] |. 
Note that for circular orbits Eq. (fSOl ) becomes 



Q 



^Si + ^L 



dE 



dL, 



(51) 



[use Eq. ( |48] | and the fact that u^^* = for circular orbits]. 
(45) This condition ensures^ that circular orbits keep circular under 



' Note in particular that the proof presented in Ref. [471, which was con- 
cerned mainly with radiation reaction, applies also to the case of the hy- 
drodynamic drag. Note also that the resonance condition which was found 
in Ref. |47] as the only possible case that could give rise to a non-circular 
evolution for an initially circular orbit is never satisfied in a Kerr space- 
time t4il. 
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the hydrodynamic drag and in the adiabatic approximation, as 
it happens for radiation reaction. 

Finally, let us note that the rates of change of E, L^ and 
Q [Eqs. (im, (|44l), ^} and ^5^] go smoothly to zero as 
the velocity of the satellite relative to the fluid goes to zero. 
This is easy to check using the fact that, when v approaches 
zero, pf^'^ = 0{v) [cf. Eqs. ^ and dSTj], p'X^"^ -^ 0, 
72_i = o{v^), Ur = 0{v), U0 = 0{v), eU-ML, = 0{v) 
and U — E — 0{v'^), and using the fact that x keeps finite 
in this limit [in particular, from Eqs. dSTT i and ( [32] i it follows 
Xt = 0{v), X4> = OiM xe = 0{l) and Xr - 0{l)l Note 
that this is indeed the result that one would expect. First of 
all, a body comoving with the fluid clearly does not experi- 
ence any dynamical friction and the only active mechanism is 
accretion. The body then accretes mass and consequently en- 
ergy and angular momentum (because the fluid carries a spe- 
cific energy and a specific angular momentum). However, the 
dimensionless constants of motion E, L^ and Q entering the 
geodesic equation cannot change because of the weak equiv- 
alence principle. Pictorially, one may think of a satellite co- 
moving with a gaseous medium. Consider a sphere centered 
in the satellite, with radius small enough to ensure that the 
gas contained in the sphere has approximately the same ve- 
locity as the satellite. Suppose now that all the gas in this 
sphere is accreted by the satellite. The velocity of the satellite 
will clearly be unaffected, because of the conservation of mo- 
mentum: for the weak equivalence principle this is enough to 
ensure that the orbit of the satellite will be unaffected, in spite 
of its increased mass. 



respect to the black hole, passing through ^inc — 90 degrees 
for polar orbits. 

In order to fix the initial conditions of a geodesic, let us first 
parametrize it with the Carter time A, which is related to the 
proper time by 1461 



— = E 
d\ 



(53) 



This is a very useful choice because it makes the geodesic 
equation separable 14611 : 



Vr{r). 



dt 
dX 



Vt{r,9), 



^1 ^vm. 






(54) 



with 



Vt{r,0)/M^ = 

(f2+a2)2 



E 



A 



-r4 



-2 • 2 

a sm 



dL^ 1 



Vr{r)/M 
E{f^ + d^) - dl. 



-A 



A 



(L, - dEf + Q 



(55) 



(56) 



V0{9)/Ar =Q-Li cot-' d - 5^(1 - E^) cos^ 0, 



B. The adiabatic approximation 

At the heart of our approach is the calculation of the 
changes of the orbital parameters experienced by Kerr 
geodesies as a result of the hydrodynamic drag, and their com- 
parison with the corresponding changes introduced by radia- 
tion reaction. To this purpose, let us recall that up to initial 
conditions Kerr geodesies can be labeled by a set of three pa- 
rameters, the semi-latus rectum p, the eccentricity e and the 
inclination angle Omc- These are just a remapping of the en- 
ergy, angular momentum and Carter constant introduced in 
Sec. HITAI and are defined as iSltj 



P = 



''a + ?■? ' 


e — 


?'a + Tp 


= i-^^ 


nin 





(52) 



where r^ and r^ are the apastron and periastron coordinate 
radii, flmin is the minimum polar angle reached during the 
orbital motion and D — \ for orbits corotating with the 
SMBH whereas Z? = — 1 for orbits counter-rotating with re- 
spect to it. Note that in the weak-field limit p and e correspond 
exactly to the semi-latus rectum and eccentricity used to de- 
scribe orbits in Newtonian gravity, and that 0-^^^ goes from 
^inc — for equatorial orbits corotating with the black hole to 
^inc = 180 degrees for equatorial orbits counter-rotating with 



(57) 



F0(r, 0)IM = Z^ csc^ 0^^dE 



f'^+d'^ ^\ a^L:, 



A 



-1 - 



A 



(58) 



This means, in particular, that the r- and 6'-motions are peri- 
odic in A. The initial conditions of the geodesic can then be 
characterized by the values to and 0o of the coordinates t and 
(f) when A = 0, the value A^o of the Carter time nearest to 
A = at which r(A,.o) = Tp, and the value Aenof the Carter 
time nearest to A = at which 6* (Ago) = ^min 14211 ■ 

Let us fix the geodesic under consideration by choosing the 
parameters p, e and 0i^c so as to obtain a bound and stable or- 
bit (see Ref. [49] for details) and by choosing the initial con- 
ditions as described above. One could in principle use Eqs. 
( |43l ), (l44l i and ( |47l ) [or ( |50] l1 to compute the rates of change 
of E, Lz and Q due to the hydrodynamic drag as a function 
the Carter time A. However, because the timescale of the or- 
bital evolution due to the interaction with the torus is much 
longer than the orbital period, we can apply the adiabatic ap- 
proximation and compute instead the averages of E, L^ and 

Q over times much longer than the orbital periods. This ap- 
proximation is routinely adopted when studying the effect of 
radiation reaction on EMRIs |42, 43l|4J,|451], and it is easy to 
implement when one considers instead the effect of the hydro- 
dynamic drag. Denoting respectively with ( )( and ( )\ the 



average over an infinite coordinate time and the average over 
an infinite Carter time, we can write Hal 



Wt 



{dt/d\)x 



(59) 



where ^ is a place-holder for either E, L^ or Q. 

Using now Eq. ^ (with the assumption that the torus 
is symmetric with respect to the equatorial plane) in Eqs. 
(l43T l, (l44l i and (|47]i [or dSOll], it is easy to show that dE/dX, 
dLz/dX and dQ/d\ depend, once fixed E, Lz and Q, only 
on the r and cos-^ 6 of the geodesic under consideration - i.e., 
r = r{X, Xro) and cos^ 9 — cos^ 9{X, Ago) - and on the sign 
of w^^*, which we will denote by e^. [The dependence on this 
sign arises because of the terms due to the radial drag, as the 
quantity w^^* appearing in the definitions of cr and x can be 
expressed in terms of r and cos^ 9 using Eqs. (l53Tl-(l58ll only 
up to such a sign.] Similarly, dt/dX is given by the geodesic 
equation (|54] | and depends, once fixed E, Lz and Q, only on 
the r and cos^ 9 of the geodesic [cf. Eq. (I55]l1. 

Using now the fact that the r- and 6'-motions are periodic 
when expressed in the Carter time, we can expand the func- 
tions dE/dX, dLz/dX and dQ/dX and dt/dX appearing in 
Eq. (|59] l in a Fourier series. Noting that the oscillating terms 
average out, one can then write these equations using only av- 
erages of these functions over the r- and 6'-periods. More pre- 
cisely, writing the r- and ^-motions as r{X, A^o) = r(A— A^o) 
and 9{X, Xeo) = 9{X — Xeo) (where we have denoted with a 
"hat" a fiducial geodesic having the same E, Lz and Q as the 
geodesic under consideration and Xro — Xeo — 0), using the 
fact that dE/dX, dLz/dX and dQ/dX depend on r, cos^ 9 and 
er, and using the fact that dt/dA depends only on r and cos^ 9, 
we can easily write Eq. ( |59] l as ||42I1 



pAe/4 



m 



/o "■ dXr Jq" dXe d*/dA(f(A^), cos2 9{Xe), e^) 



rAe/4 



J^'- dXr j;"^ dXg Vt{f {Xr),C08^ 9{Xe)) 



(60) 



where ^ is again a place-holder for either E, Lz or Q. Note 
that here A^ and Ae are the r- and 6'-periods and that dE/dX, 
dLz/dX and dQ/dX and dt/dX are expressed using Eqs. ( |43] ). 
dHi, (113 [or dSil], dH) and (|5l]l-(|58]l. Using now the defini- 
tions of cr and x [Eqs. (ISTl i and (|32]|1. it is easy to check that 
the changes of E, Lz and Q arising from the radial drag aver- 
age out in the above equation because of the presence of the 
sign er (in particular Xt,X4>jX0 oc er)- As a result, one can as- 
sume pjj^g = Oab initio when computing Eq. ( |60] l and benefit 
from another small simplification since, as we have already 
mentioned, dE/dX, dLz/dX and dQ/dX would depend only 
on r and cos^ 9 if it were not for the radial drag, which brings 
in the dependence on e^. With this assumption, all of the 
integrals appearing in (|60] | can therefore be performed over 
Xr G [0, Ar/2] rather than over Ar E [0,Ar]. Finally, note 
also that the rates of change ( |60] l do not depend on the initial 
conditions A^o and Xeo of the geodesic. 

In order to reduce Eq. ( |60b to a form suitable for numeri- 
cal integration, we can express our fiducial geodesic with the 



phase variables i/j and x, defined by 

P 



^bc 



1 + e cos tp 
cos9{x) = z- cosx . 



(61) 
(62) 



Note that ip and x change by 27r during respectively an r- and 
a 0-period. Inserting then these definitions into the geodesic 
equation (HDl one gets ll45ll49ll50ll 



di; 
dA 



P 



1-e- 



■vW) 



dx 
dX 



My/l3{z^ — z_ cos^ x) 



(63) 
(64) 



where 



I o fi ^2 l^e^^ ^1 , „„„, 




V p J 




1 — e'^ p p'^ 



{P + LI + Q) 



l3^a^{l~E^) 



z+ = 



2P 



, (65) 
(66) 

(67) 



with p ~ p/M. Note that dip/dX and dx/dX, differently from 
dr/dX and d9/dX, are non-zero at the inversion points of the 
r- and 0-motions, making ip and x very useful for numerical 
integration. 

Changing the integration variables Xr and Xg to tp and x, 
Eq. ( |60b becomes 



W2 d^/dX\p... ,{r{^P),cosH{xm'e^) 

dip I dx X 

Jo 



P^J{ip)P{z+ - z_ cos2 x) 



^^ r'^^ Vt{f{^),cosH{x)){\-e') 

Jo Py/j{lp)P{z+ -~ Z_ C0S2 x) 



(68) 



Note that the two-dimensional integrals involved in these ex- 
pressions can be easily computed numerically (e.g., iterating 
Romberg's method prj) once fixed the orbital parameters p, 
e, 9inc of the geodesic under consideration. 



IV. RESULTS 

In this section we will consider constant-^ tori around Kerr 
SMBHs and compare their influence on EMRIs with that of 
gravitational wave emission (i.e., radiation reaction) in the 
adiabatic approximation. In particular, we will compute the 
rates of change ( |68] l of the energy, angular momentum and 
Carter constant due to the hydrodynamic drag, for circular 
equatorial orbits (Sec. IIV Ab and for generic (inclined and 
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Model 


a 


A/t/M 


rout/M 


K (COS) 


£/M 


nn/M 


reenter /Af 


Pcentor (g/cm=^) 


poavg (g/cm=') 


Al 


0.900 


0.100 


1.000x10^ 


4.198x10^^ 


2.6324500536 


1.73246 


3.60963 


4.060x10"=^ 


1.475x10"" 


A2 


-0.900 


0.100 


1.000x10^ 


4.189x10^^ 


4.7567317819 


5.65700 


15.58890 


3.992x10"^ 


1.476x10"" 


Bl 


0.998 


0.100 


1.000x10^ 


4.200x10^^ 


2.0894422310 


1.09144 


1.56484 


1.868x10"'' 


1.474x10"" 


B2 


0.500 


0.100 


1.000x10^ 


4.195x10^2 


3.4141929560 


2.91425 


7.16458 


1.331x10-"^ 


1.475x10"" 


B3 


0.000 


0.100 


1.000x10^ 


4.192x10^2 


3.9999599993 


4.00008 


10.47174 


7.355x10"'' 


1.475x10"" 


B4 


-0.500 


0.100 


1.000x10^ 


4.190x10^^ 


4.4494291313 


4.94962 


13.39547 


5.034x10"'' 


1.475x10"" 


B5 


-0.998 


0.100 


1.000x10^ 


4.188x10^2 


4.8269302324 


5.82521 


16.11218 


3.796x10"'* 


1.476x10"" 


CI 


0.900 


0.100 


1.000x10=* 


3.997x102" 


2.6319080229 


1.73318 


3.60622 


4.323xl0"2 


1.507x10"='' 


C2 


-0.900 


0.100 


1.000x10^ 


3.607x102" 


4.7490561067 


5.67540 


15.49760 


4.775x10"=' 


1.573x10"^ 


C3 


0.900 


0.100 


l.OOOxlO'' 


4.170x1021 


2.6324007478 


1.73253 


3.60932 


1.295x10-=* 


1.479x10"=* 


C4 


-0.900 


0.100 


LOOOxlO-* 


4.103x1021 


4.7560304461 


5.65866 


15.58057 


1.297x10"'' 


1.489x10"=* 


C5 


0.900 


0.100 


1.000x10*^ 


4.201x102=* 


2.6324549842 


1.73246 


3.60967 


1.283x10"'* 


1.475x10"" 


C6 


-0.900 


0.100 


l.OOOxlO'' 


4.199x10^=* 


4.7568019526 


5.65683 


15.58974 


1.258x10"'^ 


1.475x10"" 



Table I: Models analyzed in this paper: all of them have M — 10^ Mq, F — 5/3 and are filling exactly their outermost closed equipotential 
surface (i.e., they have AW = 0). All the parameters are defined in Sec. [ill except the average rest-mass density poavg ~ Mto/V, where 
V — f „ y/—gd X. Note that the specific angular momentum of the torus needs to be tuned with high accuracy in order to obtain large 
outer radii such as those considered in these models, and for this reason we report l/M with 10 decimal digits. 



eccentric) orbits (Sec. lIVBI l. Since E, L^ and Q can be ex- 
pressed analytically as functions of the orbital parameters p, 
e and 0inc [49], it is then easy to compute the rates of change 
dp/dt, de/dt and d9ij-^c/dt due to the hydrodynamic drag. For 
the same orbits, we will consider also the radiation reaction, 
for which we will use the kludge fluxes dE/dt, dL^/dt and 
dQ/dt of Ref. |52] to compute dp/dt, de/dt and d9inc/dt. 
Note that these kludge fluxes are a good approximation to 
the fluxes computed rigorously in the adiabatic approximation 
with the Teukolsky-Sasaki-Nakamura formalism ll45ll53ll . In 
fact, since they are based on a Post-Newtonian expansion cor- 
rected with fits to fluxes computed with the Teukolsky-Sasaki- 
Nakamura formalism for circular orbits, these kludge fluxes 
are accurate within 3% for circular orbits and their accuracy 
is expected to be within 10 — 15% also for generic orbits with 
p > 6AI. Moreover, they are expected to be off at most by 
25 — 30% even for smaller values of the semi-latus rectum p 
(cf. Ref. 152], Table I). 

The mass of the SMBH is fixed to M ^ IO^Mq while its 
spin parameter a ranges from — 0.998A/ to 0.998M (note that 
\a\ = 0.998A/ is a reasonable upper limit for the spin attain- 
able as the result of mass accretion [54] or binary black-hole 
mergers 11551 Dal ), and the mass of the satellite black hole is 
instead m = lAf©. The constant-^ torus is assumed to be 
composed of an isentropic monatomic gas (i.e., a F = 5/3 
polytrope) and is considered to be exactly filling its outermost 
closed equipotential surface ( A14^ — 0), so as to present a zero 



accretion rate Af onto the SMBH.^ Once assumed /SW = 
and F = 5/3, the specific angular momentum of the torus 
is uniquely fixed by choosing the outer radius. A reason- 
able outer radius for a realistic accretion disk is given by 
''out = 10^ Af [141 . and this is indeed the value that we will 
use for most of our analysis, although we will briefly consider 
also different values for rout in order to study the impact of 
this parameter on the final results. The poly tropic constant 
K of the equation of state is finally fixed by the requirement 
that Aft = O.lAf . While this could be a reasonable value for 
the mass of a realistic accretion disk in AGNs OTJIl . we will 
see that our results scale proportionally to Aft, thus allowing 
one to extrapolate them easily to the case A^t = Af , which is 
certainly an astrophysically plausible value, but one for which 
our test-fluid approach is no longer valid. 

We should stress, however, that our results, when expressed 
in terms of the dimensionless orbital parameters p/M, e and 



' While realistic thick disks are generally expected to accrete onto the 
SMBH, these configurations are clearly non-stationary and cannot there- 
fore be reproduced within our framework. However, it is easy to show that 
if one cuts-off a torus solution with M > at r = reusp, the effect of 
the satellite-torus interaction will be enhanced with respect to the M = 
solution having the same mass and outer radius: the choice Af = is thus 
useful to obtain at least a lower limit for the effect of the satellite-torus 
interaction on EMRIs. 
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Figure 1: The top plot shows the absolute value of the ratio f3r = {dr/dt)hydTo/(d.r/dt)Gw between the rates of change of the orbital radius 
due to hydrodynamic drag and radiation reaction, for circular bound stable orbits in the equatorial plane, as a function of the radius r. For 
graphical reasons r is considered positive for orbits in the positive (ji-direction ("prograde orbits" i.e., corotating with the torus) and negative 
for those in the negative ^-direction ("retrograde orbits" i.e., counter-rotating with respect to the torus). Note that all the curves of this figure 
(including those of the middle and bottom plots) are terminated at the (prograde or retrograde) ISCO. The middle plot compares the velocity 
V of the satellite in the rest frame of the fluid with the sound velocity Vs, while the bottom plot shows the energy density of the torus. The 
curves refer to the models Al and A2 of Table|l] which are labeled here with the spin parameter a of the SMBH. Note that in all the plots the 
vertical axis is drawn in logarithmic scale. As such, the vertical asymptotes appearing in these plots actually correspond to a zero value for the 
quantity under consideration. 



^inc, are approximately independent of the mass M of the 
SMBH and of the mass m of the satellite black hole (provided 
that Mt/M and rout/M are maintained constant). Indeed, 
since the ratios between the rates of change dp/dt, de/dt and 
dOinc/dt due to the hydrodynamic drag and radiation reaction 
are of course dimensionless, it is not restrictive to fix M = 1, 
because this simply corresponds to choosing a system of units. 
Note in particular that this means that systems with different 
m and M but equal mass ratio ra/M give exactly the same 
ratios between the rates of change dp/dt, de/dt and d9inc/dt 
due to the hydrodynamic drag and radiation reaction. More- 
over, these rates of change are proportional to m (in the case 
of the hydrodynamic drag this can be seen from Eqs. dZTJ l. 
(l33T l. (|3J] |. (|43] |. (Ell and (|47]i, while in the case of radiation 



reaction see for instance Ref. 115 2|] ). so this dependence on m 
cancels out when taking the ratio. The only dependence on 
m arises from the cutoff r^in ^ 2m (1 + w^)/w^ appearing 



in Eq. (|35] |, but this dependence clearly comes about only for 
supersonic velocities and is a logarithmic one. As such, the 
results which we present in this paper, although derived in 
the case of m = lAf© and M — 10^ Af©, are also valid for 
m = O.lAfo and M = 10^ Mq (exactly) or for m = lAf© 
and M — 10^ Af© (exactly for subsonic motion, and approxi- 
mately - with an error comparable with those affecting the fit 
(|35] l or the kludge fluxes - for supersonic motion). 

In all of our analysis we will focus on the region close to 
the SMBH (r < 50A/), which contains only a small fraction 
of the mass of the torus (e.g., in the case of the model Al of 
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Table m the mass contained in a radius r — 50 Af amounts to 
about 2.9 X 10^^ Mt, and this fraction scales approximately as 

— 3/2 

r^^^ when considering tori with different outer radii). This 
is the region relevant for gravitational-wave experiments like 
LISA. In particular, an EMRI's signal is expected to be de- 
tectable by LISA when its frequency (which is twice the or- 
bital frequency) increases above w 2 mHz (below this fre- 
quency, in fact, there is a strong unresolvable foreground noise 
due to double white-dwarf binaries in our Galaxy ISTfl). This 
translates into a distance from the SMBH of r « lOM for 
M = IO^A/q, and to r w 45Af for M == IO^A/q. 

It should be noted, however, that the amplitude of an 
EMRI's signal scales with the distance from the SMBH: for 
a circular orbit of radius r, the Keplerian frequency is 27ti' = 
Af i/V(r3/2 ± aA/i/2) ^ ^1/2/^3/2 ^j^^j ^^^ ampHtude of 

the signal is /i ~ (m/L>)(27r:/Ai")2/3 ~ {m/D){M/r) 1,58,1 , 
where D is the distance from the observer to the source. 
As such, an EMRI ai'ound a IO^Mq SMBH will have a 
gravitational-wave amplitude that at r ~ 45 Af is about 10 
times smaller than at r ~ 5A/. Therefore, to see the de- 
tails of the waveforms at r ^ 45A/ the source must be 
~ 10 times closer to us, which translates into a detection vol- 
ume decreased by a factor ~ 1000. Nevertheless, this de- 
crease of the detection volume may be compensated (at least 
partly) by the fact that the event-rate estimates consider only 
EMRJs in the strong-field region of the SMBH, even when 
M — 10^ A/0 [3]. As such, since EMRIs in the early part of 
the inspiral are more numerous than those in the strong-field 
region, one expects to see a number of these events larger 
than the naive estimate given by the rate expected for strong- 
field EMRIs around a 10^ A/© SMBH divided by the detection 
volume decrease factor ^ 1000. Of course, the event rates 
could be even larger if the satellite were a black hole with 
m ^ IOOATq, because the amplitude of the signal is propor- 
tional to m, but too little is presently known about these ob- 
jects to draw any sound conclusions (see for instance Ref. ||59I1 
for a review on intermediate-mass black holes as possible 
sources for LISA). 



A. Circular equatorial orbits 

The evolution of circular equatorial orbits is very simple 
in the adiabatic approximation. As mentioned in Sec. IIII Al 
both the radiation reaction and the hydrodynamic drag main- 
tain circular orbits circular and, due to the symmetry of the 
Kerr spacetime and of the torus with respect to the equato- 
rial plane, equatorial orbits will remain equatorial. There- 
fore, the evolution of circular equatorial orbits under both 
radiation reaction and hydrodynamic drag can be character- 
ized with only one quantity (the rate of change of the radius 
dr/dt), to which the rates of change of the energy and an- 
gular momentum, dE/dt — (dE/dr) (dr/dt) and dL^/dt — 
[dLz/dr) (dr/dt), are proportional. {dQ/dt is instead identi- 
cally zero for equatorial orbits.) Moreover, one does not need 
to compute the infinite-time averages ( |68] l, because the rates of 
change of E, L^ and Q due to the hydrodynamic drag, given 
by Eqs. (|43]l, ^^ and ^} [or ^5^], are akeady functions of 



the orbital radius alone. [Note also that the 4-vector x reduces 

to X = ~dr/y^.] 

The ratio between the rates of change of the orbital ra- 
dius due to hydrodynamic drag and radiation reaction is 
a convenient measure of the "efficiency" of the hydro- 
dynamic drag. Defining this quantity simply as (3f = 
(dr/dt)\^yi^o/(dr/dt)Q\f^, we show in the top plot of Fig. [T] 
the absolute value of (3^ as a function of the radius r of cir- 
cular equatorial bound stable orbits. The two curves refer to 
models Al and A2 of Table U and are labeled with the spin 
parameter a of the SMBH. Note that in order to present all the 
data in a single plot, a positive r refers to orbits rotating in the 
positive (/)-direction ("prograde orbits" i.e., corotating with the 
torus), while a negative r refers to orbits rotating in the nega- 
tive (/)-direction ("retrograde orbits" i.e., counter-rotating with 
respect to the torus). The middle plot compares the velocity v 
of the satellite in the rest frame of the fluid with the sound ve- 
locity Vs, while the bottom plot shows instead the energy den- 
sity p of the torus. Note that in all the plots the vertical axis is 
drawn in logarithmic scale. As such, the vertical asymptotes 
appearing in Fig. [T](as well as in Figs. |2]and|3] which are in 
logarithmic scale too) actually correspond to a zero value for 
the quantity under consideration. 

Note that if the torus is corotating with the black hole, the 
radius of the innermost stable circular orbit (ISCO) is always 
larger than the inner radius of the torus, both for prograde and 
retrograde orbits. For prograde orbits, this immediately fol- 
lows from the considerations of Sec. |II] (since our tori have 
AVK = and £ms < ^ < -^mb, we have rin = rcusp < ^ns, 
and Tins is exactly the radius of the prograde ISCO), while 
for retrograde orbits it is sufficient to note that the retrograde 
ISCO is located at a larger radius than the prograde one. Bear- 
ing this in mind, it is then easy to understand why none of the 
quantities plotted in Fig. [T]for model Al (a = 0.9A/) goes to 
zero when approaching the SMBH: although the density, the 
velocity of the satellite relative to the torus, the sound veloc- 
ity and (dr/di) hydro are exactly zero at the inner radius of the 
torus, this radius is smaller than that of the ISCO and therefore 
no bound stable orbits exist there. 

If instead the torus is counter-rotating with respect to the 
black hole {i.e., a < 0), the radius of the ISCO is larger than 
Tin for prograde orbits (this follows again from rjn = rcusp < 
Tins), but it is not possible to conclude that the radius of the 
ISCO is larger than r^n also for retrograde orbits. In fact, the 
ISCO counter-rotating with respect to the torus (i.e., the "ret- 
rograde" ISCO) is corotating with the black hole and thus lies 
at a radius smaller than the "prograde" ISCO. Indeed, for the 
model A2 (a = — 0.9Af) considered in Fig. [T] the retrograde 
ISCO is at a radius smaller than r^n- As a consequence, the 
density, the sound velocity and ((ir-/(it)hvdro for model A2 go 
to zero when the radius of the retrograde orbits decreases, be- 
ing in fact zero at the inner edge of the torus. (Of course, the 
velocity of the satellite relative to the fluid does not go to zero 
when the radius of the retrograde orbits decreases, because the 
satellite and the torus are rotating in opposite directions.) 

As it can be seen in the top plot, the ratio \Pf\ is larger for 
the retrograde orbits than for the prograde ones. The reason 
can be easily understood from the middle plot, which shows 
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Figure 2: The same as the top plot of Fig. [T] but for the models B1-B5 of Table |l] which are labeled here with the spin parameter a of the 
SMBH. 
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Figure 3: The same as the top plot of Fig. [T] but for the models A1-A2 and C1-C6 of TableU which are labeled here with the spin parameter a 
of the SMBH and the outer radius rout of the torus. 
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that the retrograde motion is always supersonic. As such, the 
long-range drag, which increases significantly when passing 
from the subsonic to the supersonic regime [cf. Eq. (|35]|1, en- 
hances the torus-satellite interaction for the retrograde orbits. 
From the middle plot one can also note that relativistic veloci- 
ties {v > 0.6) are reached in the case of retrograde orbits very 
close to the SMBH, thus further enhancing the hydrodynamic 
drag because of the relativistic correction factor 7^(1 + w^)^ 
appearing in Eq. ( l33T l. However, we should note that when v 
reaches its maximum value {i.e., v ~ 0.8 for model A2) the 
effect of the relativistic correction factor on the hydrodynamic 
drag is hindered by the small value of the density, which goes 
to zero at the inner edge of the torus {cf. the bottom plot). 

As it can be easily understood from the formulas reviewed 
in Sec lUl a change in the polytropic constant k leaves all the 
parameters of the torus unchanged, except the energy density 
p, the rest mass density po, the pressure p [all of which scale 
proportionally to k^^^^^^^^] and the total mass-energy and 
rest-mass, Alt and Afto (which scale proportionally to p). For 
this reason, the bottom plot of Fig. [T] scales linearly with the 
mass of the torus {cf. the label of the vertical axis). Noting, 
from the formulas of Sec. IIII Al that the rates of change of the 
E, Lz and Q are proportional to the energy density p, the same 
scaling applies to the top plot. This is a very useful feature, 
because although the value used for the figures of this paper 
- i.e., Mt = O.IM - could be a plausible mass for the torus, 
very little is known about these objects and larger or smaller 
masses may be possible. In general, a different mass Alt could 
have important effects. For instance, extrapolating to the case 
Alt — M, in which our test-fluid approximation is no longer 
valid, the ratio |/3j. | would be enhanced by a factor 10, and for 
orbits counter-rotating with respect to the torus with r > 40M 
the effects of hydrodynamic drag and radiation reaction would 
become comparable.^ 

It is also worth pointing out that at the center of the torus 
(dj'/di) hydro changes sign forprograde orbits, being negative 
forr > reenter and positive for r < reenter- [{dr/dt)Gw is in- 
stead always negative because gravitational waves carry a pos- 
itive amount of energy away from the source.] This change of 
sign corresponds, in the top plot of Fig. [T] to the zero value for 
\f3r\. This behavior comes about because, although the density 
reaches its maximum at the center, the motion of the fluid is 
exactly Keplerian (geodesic) there, and the relative velocity 
of the satellite is therefore exactly zero {cf the middle plot). 
This means in particular that E = — Mj^' = — uj"*""* = U, 

which together with Eq. ( l43T l and xt — gives E — and 
therefore f — for prograde orbits at r = reenter- Moreover, 
if r > reenter the spccific angular momentum of the satellite is 
larger than that of the torus {cf for instance Ref. II33[1 . Fig. 5), 



' Note, however, that even values of \f3r\ less than 1 can produce fea- 
tures detectable by LISA, because the dephasing time scales as ~ 
I'^Gwl"^ fedl - For instance, if |/3y,| Ri 0.1 the dephasing time between 
a waveform with only radiation reaction included and one with also the ef- 
fect of the hydrodynamic drag is expected to be only ~ 3 times larger than 
the dephasing time between waveforms with and without radiation reaction 
included. 



and therefore the satellite is slowed down by the interaction 
with the fluid {i.e., f < 0). On the other hand, if r < reenter 
the specific angular momentum of the satellite is smaller than 
that of the torus, and the satellite is speeded-up {i.e., f > 0). 

Fig. |2] shows the absolute value of (3f as a function of 
the radius r of circular equatorial bound stable orbits for the 
models B1-B5 of Table U whose spin parameter a goes from 
— 0.998Af to 0.998A/. As it can be seen, the situation is qual- 
itatively very similar to the one presented in the top plot of 
Fig. [T] In particular, the effect of the hydrodynamic drag can 
be comparable to that of radiation reaction, but only if we ex- 
trapolate to Mt = M and, even in that case, only for orbits 
counter-rotating with respect to the torus and with r > 40Af . 

We can also note that the effect of the spin a on the results 
is negligible, except for the prograde orbits between the cen- 
ter and the ISCO, for which \/3r\ decreases as a increases. The 
reason for this can be easily understood by considering a satel- 
lite moving on a prograde circular equatorial orbit between the 
center and the inner edge of the torus, and by recalling that the 
difference between ^mb and £nis represents an upper limit for 
the deviation of the specific angular momentum of the satellite 
away from that of the torus (see Ref. ||33I1 . Fig. 5). Because 
this deviation regulates the exchange of angular momentum 
between the torus and the satellite [cf. Eq. ( l44l l. where Xct> ^ 
for circular orbits] and thus the rate of change of the orbital ra- 
dius, f3r must go to zero if ^,„b — ^ms goes to zero. Since it is 
easy to verify that ^nib — ^ms ^ as a ^ AI [cf. Eq. jTSJ], it 
is natural to find that \f]r\ decreases as a increases. 

Finally, in Fig. [3] we plot again \(3f\, but for models Al- 
A2 and C 1 -C6 of Table U in which we have considered dif- 
ferent values for the outer radius rout of the torus, ranging 
from lO^Af to 10^ Af. The reason for this is that although 
''out ^ 10^ Af is a plausible value for the outer radius, little 
is known about the size of astrophysical accretion disks and 
larger or smaller outer radii may also be possible. As it can 
be seen from Fig. |3] a different outer radius will have signifi- 
cant effects for prograde orbits with r > 20A/ and retrograde 
orbits with r > lOA/, for which the effect of the hydrody- 
namic drag can become comparable to that of radiation re- 
action. In general, \f3r\ progressively increases as the outer 
radius is decreased. This is rather simple to explain: decreas- 
ing rout while keeping Alt constant amounts to increasing the 
average rest-mass density and hence the hydrodynamic drag. 
When considered from this point of view, the uncertainty on 
the value of rout has an effect opposite to the the uncertainty 
about the mass of the torus: a decrease of rout (or an increase 
of Alt) induces an increase of |/3r|. For circular orbits, this 
overall uncertainty can be easily modelled in terms of a sim- 
ple scaling of the type 



'lO^Ai'X^''^ 

/3r(?'out) ~ (3f\5 [ — 1 



rout 



(69) 



where /3r|5 is the efficiency for rout — 10 AI. Note that the 
scaling power is not 3 as one may naively expect. This is 
because (3f. is most sensitive to the changes of the rest-mass 
density in the inner part of the torus and this does not scale 
simply as roui 
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Figure 4: log^Q |/3p | for inclined orbits witii e = 0. 1 as a function of the semi-latus rectum p and of the inclination angle 6inc- The figure refers 
to the model Al of Table|l] but a larger (smaller) mass A/t for the torus would simply increase (decrease) the absolute ratio |/3p| by a factor 
Mt/(0.1Af). The dashed line marks the edge of the torus. 



B. Generic orbits 



We will now extend the analysis of Sec. IIVAI to bound 
stable generic (inclined and eccentric) orbits. Although such 
an extension is in principle straightforward using the formu- 
las introduced in Sec. 1111 B I the space of parameters and re- 
sults which one needs to examine greatly enlarges. Not only 
are generic orbits characterized by the three parameters p, e, 
^inc defined by Eq. ( |52] |. but one also needs to consider three 
quantities describing the evolution of each single orbit in the 
parameter space i.e., the rates of change dp/dt, de/dt and 

ddinc/dt. 

To simplify our analysis, we will focus mainly on model Al 
of Table U which could be a representative example of an as- 
trophysical torus in an AGN, and then examine how the rates 
of change dp/dt, de/dt and dOmc/dt due to the hydrodynamic 
drag compare to those due to radiation reaction throughout 
the space of parameters (p, e, 6'inc)- The considerations that 
we will draw for model Al can, however, be extended sim- 
ply to the cases of different masses and radii for the torus. As 
in the case of circular orbits, in fact, a larger (smaller) mass 
Mt for the torus when Tout is held constant would simply in- 
crease (decrease) the rates dp/dt, de/dt and ddinc/dt due to 
the hydrodynamic drag by a factor A/t/(0.1A/). This scaling 
is exact (as long as the torus is not self-gravitating) and comes 
about because the rates of change of E, Lz and Q (and conse- 
quently those of p, e and 6'iiic) are proportional to the energy 
density p ex Mt- Similarly, variations of rout will result in an 
effect which is similar to the one discussed for Fig. [3] in the 



case of circular orbits [cf., eq. (|69]l1, as we will see at the end 
of this section. 

All of the results presented in this Section have been com- 
puted by integrating numerically Eq. ( |68] | using an iterated 
Romberg method |51], with a typical accuracy, depending on 
the parameters of the orbit under consideration, of 10^^-10^'* 
and never worse than 4x10" 



-3 10 



We start by analyzing in detail orbits with eccentricity 
e = 0.1, and we will then study the effect of a different 
eccentricity on the results. In Figs. |4]|5]|6]the color-code and 
the contour levels show the base- 10 logarithm of the "efficien- 
cies" of p, e and 6'inc, i-e., Pp = |(dp/di)hydro/(c^p/di)Gw|, 
/3e = |(de/(ii)hydro/(rfe/di)Gw| and ^g._^^ = 
|(c^6'inc/rfi)hydro/(rf^inc/rfi)Gw|, as functions of p and 
^inc- The dashed line marks the "edge of the torus" i.e., 
the location in the {p, 6'inc)-plane of the orbits having 
^inc = ^inc,t(p), where 9i-ac,t{r) is the function giving the 
angle between the surface of the torus and the equatorial 
plane in terms of the radius r. The blank part on the left 
portion of these figures refers to the region where no bound 
stable orbits exist, and we will refer to the line marking the 
boundary of this region as the separatrix Ii61i] 



'" Note that the accuracy of the numerical integration is certainly adequate, 
because it is considerably better than the errors affecting the fit (35) as 
well as those affecting the kludge fluxes that we use to study the effects of 
radiation reaction. 
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Figure 5: The same as in Fig. |4l but for logj^Q |/?e 



Each figure has been obtained by computing the quantity 
under consideration using Eq. (l68T l for ^^ 5 x 10^ orbits ir- 
regularly distributed in the {p, Sine) -plane, and then linearly 
interpolating on a grid of 1500 x 1500 nodes using a Delau- 
nay triangle-based method. The gridded data obtained in this 
way has been used to draw the contour levels. Not surpris- 
ingly. Figs, m |5] m show somewhat the same trend as the 
results presented in Sec. IIV Al for circular equatorial orbits, 
with the effect of the torus becoming comparable to that of the 
radiation-reaction far away from the black hole and becoming 
instead negligible in the strong-field region of the black hole. 
However, these figures present also a variety of features that 
we will now analyze in detail. 

Fig. m for instance, shows \(3p\ and indicates that the ef- 
fect of the hydrodynamic drag is larger for orbits with high 
inclination, for any given semi-latus rectum p. This is simply 
due to the fact that orbits with 9inc > 90 degrees are retro- 
grade with respect to the torus, and the velocity of the satel- 
lite relative to the fluid can easily become supersonic. Indeed 
this effect is visible also in the figures of Sec. IIV Al (We 
recall that in those figures the retrograde orbits are mapped 
to negative values of the orbital radius r.) The transition be- 
tween the subsonic and the supersonic regime is marked by 
the sharp bend of the contour levels of Fig. |4]at ^jnc ~ 40 
degrees. This bend corresponds to the passage from orbits 
which are always subsonic (the orbits with ^jnc smaller than 
the inclination angle at which the bend is located) to orbits 
which are supersonic at least for a part of their trajectory (the 
orbits with 9inc larger than the inclination angle at which the 
bend is located). Another small dip is hardly noticeable in the 
contour levels at inclination angles 9inc just smaller than the 



edge of the torus (and smaller than 90 degrees, corresponding 
therefore to prograde orbits); this feature corresponds to the 
transition from orbits which are partly subsonic and partly su- 
personic ("below" the dip), to orbits which keep always super- 
sonic ("above" the dip)." From Fig. |4]one can also note that 
|/3p| becomes lower than 10^^ in anarrow "strip" atp/M < 5. 
Indeed, ((ip/di)hydro changes sign inside this "strip", being 
positive inside the region between the "strip" and the separa- 
trix and negative outside, while {dp/dt)cw is always nega- 
tive. This behavior generahzes that of circular equatorial or- 
bits, for which {dr / dt)i^ydio changes sign at the center of the 
torus {cf. Sec. IIV Al l. Also in this case, however, the very 
small values of |/3p| cannot produce an observable imprint on 
the waveforms. 

In a similar way. Fig. |5]shows the behavior of | /3e | . As it can 
be seen, the influence of the torus is again larger at high incli- 
nations than at low ones, for any fixed semi-latus rectum. Also 
in this case, this happens because the orbits counter-rotating 
with respect to the torus can easily become supersonic. As in 



" We note that in order to better understand the fine features in the con- 
tour levels, we have built an auxiliary code computing the quantities 
(dp/dt)hydro. (rfe/df)hydro and (c(cos6'inc/dt)hydro by direct integra- 
tion of Eqs. J43K I44t . j47t and j5Qt along numerically solved geodesies, 
averaging over a reasonably large number of revolutions (~ 30) for each 
geodesic. This has not only validated the results which have been used to 
build the figures and which have been obtained using Eqs. (68), but has 
also allowed us to examine in detail the behavior of the geodesies in the 
various regions of the parameter space (p, e, Sine), thus helping to inter- 
pret the complicated features of the figures shown in this paper. 
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Figure 6: The same as in Fig. [4] but for logj^Q \Pg 



Fig. m we can note the presence of a sharp bend in the con- 
tour levels at Oi^c ^ 40 degrees, due to the transition from 
orbits always subsonic to orbits partly supersonic, and a dip 
in the contour levels near the edge of the torus (at inclina- 
tions Oiric < 90 degrees), more pronounced than in Fig. H] 
but again due to the transition from orbits only partly sub- 
sonic to orbits always supersonic. Moreover, one can note the 
presence of three "valleys" where the efficiency |/3e| becomes 
very small. One ("valley 1") starts at 6'i,ic ~ 15 degrees, very 
close to the separatrix, and extends as far as the right edge of 
the figure {p/M = 50, Sine ~ 65 degrees) and beyond. A 
second valley ("valley 2") starts at the same point as valley 
1, but extends only until p/M « 12 and 6'inc ~ 60 degrees, 
where it terminates together with a third valley ("valley 3") 
starting on the separatrix at 6'inc ~ 30 degrees. Across these 
valleys, the quantity (de/(it)hydro becomes zero and changes 
sign, being negative under valley 1 and in the region between 
the separatrix and valleys 2 and 3, and positive in the rest of 
the {p, 6'inc) -plane. Conversely, the rate of change of the ec- 
centricity due to radiation reaction is always ne gati ve, with the 
exception of orbits very close to the separatrix 115211 ; this is ap- 
parent also in Fig. |5] where the narrow "strip" corresponding 
to a ratio |/3e| > 10^'^ and running close and almost paral- 
lel to the separatrix is due to a change in sign of {de/dt)^^. 
Despite this markedly different behavior of (de/di) hydro and 
{de/dt)Q^fij, Fig. |5]shows that the effect of the hydrodynamic 
drag is always much smaller than radiation reaction unless the 
semi-latus rectum of the orbit is increased to p/M > 50, 
or the mass of the torus is increased at least by a factor 10 



thus extrapolating to Mt = A/'^. However, while a larger 
semi-latus rectum increases the efficiency /3e, it also reduces 
the frequency and amplitude of the gravitational-wave signal, 
moving it to a region of low sensitivity for LISA. 

The effect of the hydrodynamic drag is somewhat stronger 
when considered in terms of the efficiency |/3g |, as it is 
shown in Fig. |6] While the qualitative behavior is similar to 
the one discussed for the two preceding figures, it should be 
noted that the effect of the torus is comparable to that of radia- 
tion reaction akeady for p/M w 32 if ^inc ~ 60 degrees, and 
the two effects remain comparable down to p/M « 20 if the 
mass of the torus is increased by a factor 10 thus extrapolat- 
ing to A/t — M. Moreover, while radiation reaction produces 
an increase in the inclination 6'inc irrespective of the orbital 
parameters pa], (d6'inc/o'0 hydro is always negative and thus 
a measurement of the evolution of the inclination angle 6'inc 
in the early stages of an EMRl could give important informa- 
tion on the presence of a torus. If such a presence were to 
be detected, it would not prevent high-precision tests of the 
Kerr nature of the SMBH being performed in the strong-field 
region, where the hydrodynamic drag becomes negligible. 

We should also note that the decrease of 9inc due to the hy- 
drodynamic drag is not surprising for orbits with 6'inc > 90 
degrees (i.e., orbits rotating in the opposite 0-direction with 
respect to the fluid), because the hydrodynamic drag clearly 
induces the orbits to rotate in the same (7!)-direction as the 



■ The test-fluid approximation of course breaks down in this hmit. 
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Figure 7: logj^Q |/3g. j shown as a function of the semi-latus rectum p and of the inclination angle ^inc for inclined orbits with selected values 
of the eccentricity, i.e., e = (solid line), e = 0.4 (dot-dashed line), e = 0.8 (dotted line). The figure refers to the model Al of TableU] but a 
larger (smaller) mass Aft for the torus would simply increase (decrease) |/3g. j by a factor Aft/(0.1Af ). The dashed line marks the egde of 
the torus. 



torus. For orbits with ^inc < 90 degrees, instead, the de- 
crease of ^inc comes directly from Eq. ( |45l ) (we recall that the 
effect of the radial drag averages out when adopting the adi- 
abatic approximation), thus following directly from the axis- 
and plane-symmetry of the system and being independent of 
the use of constant specific angular momentum tori such as 
the ones considered in this paper. Indeed, since the fluid of 
the torus does not move in the 6'-direction, Eq. ( l45b states 
that accretion conserves the momentum of the satellite in the 
6'-direction, but it also increases its mass, thus reducing the 
velocity in the 6'-direction. In addition, the dynamical fric- 
tion will contribute to this reduction by damping further the 
oscillations around the equatorial plane. 

To illustrate how the above results depend on the eccentric- 
ity, we show in Fig.|2]the efficiency |/?g | for the model Al of 
TableU but for different values of the eccentricity i.e., e = 0, 
0.4, and 0.8. (Equivalent figures could be made also for |/3p| 
and |/3e|, but we omit them here because they would be qual- 
itatively similar to Fig. |7]l. For each value of the eccentricity, 
we have computed |/3g | for ^ 4 x 10** orbits, and using the 
same technique employed for Figs. |4]|5]|6]we have drawn the 
contours corresponding to values of 0.01, 0.1, 1 and 10. Also 
in this case, a larger (smaller) mass Aft for the torus would 
simply increase (decrease) these absolute ratios by a factor 
Aft/(0.1A/). 

Clearly, many of the features in this plot have been dis- 
cussed also for the previous figures. For instance, the con- 
tour levels present sharp bends at low inclinations {i.e., B\^c ~ 



10 — 40 degrees) for e = and e = 0.4, due the transition 
from subsonic to partly supersonic orbits, whereas the transi- 
tion from partly supersonic to fully supersonic orbits causes 
the appearance of a pronounced "kink" in the e = 0.4 contour 
levels, for inclinations 6'iiic < 90 degrees just above the edge 
of the torus. The e = 0.8 contour levels, on the other hand, 
are rather smooth and less affected by the complex changes 
of regimes as the satellite interacts with the torus. Most im- 
portantly, however. Fig. [T] suggests that the conclusions drawn 
when discussing Fig. |6]for orbits with e = 0.1 are not al- 
tered significantly by a change in the eccentricity. Indeed, 
even for large eccentricities the influence of the torus on the 
evolution of 6'inc can be comparable to that of radiation reac- 
tion for p/M as small as 35 — 38, while the two effects are 
still equal at p/A/ w 23 — 24 if the mass of the torus is in- 
creased by a factor 10 thus extrapolating to A/t = M. As a 
result, a measurement of the evolution of the inclination an- 
gle 6'inc even for generic eccentric orbits could give important 
information on the presence of a torus around the SMBH. 

This conclusion is finally confirmed by Fig. [8] in which 
we show how |/3g | changes if one considers different values 
of the outer radius of the torus while keeping its mass fixed. 
More specifically. Fig. [8] shows the location in the (p, 6'inc)- 
plane of the circular orbits for which |/3g | = 1- Different 
curves refer to different values of the outer radius, and in par- 
ticular to the models Al, CI, C3 and C5 of TableU As it is 
probably obvious by now, a different mass Aft for the torus 
would simply make the curves of this figure correspond to 
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\/3g. I — Mt/{0.1M). As expected from the results of Sec. 
IIV A| a modest variation of the outer radius can easily cause 
the decrease of the inclination angle due to the hydrodynamic 
drag to be dominant over the increase due to radiation reaction 
for orbits with p/M ^ 20 or smaller. 



V. CONCLUSIONS 

SMBHs are expected to be surrounded by matter, either in 
the form of stellar disks, as in the case of "normal" galactic 
centers, or in the form of accretion disks of gas and dust, as 
in the case of AGNs. In order to assess under what conditions 
and to what extent the interaction with matter could modify 
the gravitational-wave signal from EMRls in AGNs, we have 
studied EMRls in spacetimes containing a SMBH suiTounded 
by a non self-gravitating torus. For simplicity, and in order 
to handle the equilibrium solution analytically, we have con- 
sidered a torus with a constant distribution of specific angular 
momentum, using as reference dimensions and masses those 
for the accretion disks expected in AGNs, but bearing in mind 
that these also come with rather large uncertainties. We have 
extrapolated our results also to cases in which the mass of the 
torus is comparable with that of the SMBH, although we stress 
that in this limit our test-fluid approximation for the torus is 
no longer valid. 

Overall, we have found that the effect of the hydrodynamic 
drag exerted by the torus on the satellite black hole can have 
important effects sufficiently far from the central object, and 
that these effects are qualitatively different from those of ra- 
diation reaction. In particular, if the torus is corotating with 



the SMBH, the hydrodynamic drag always decreases the in- 
clination of the orbits with respect to the equatorial plane (i.e., 
orbits evolve towards the equatorial prograde configuration), 
whereas radiation reaction always increases the inclination 
(i.e., orbits evolve towards the equatorial retrograde configura- 
tion). In the case of a system composed of a SMBH with mass 
M — 10® Mq and a corotating torus with mass Aft < M, the 
effect of the torus will be marginally observable by LISA only 
if the radius of the torus is as small as rout ~ 10"^ — 10'* A/. 
However, if the SMBH has a lower mass, EMRls will be de- 
tectable by LISA at larger distances from the SMBH, and the 
effects of a torus will be more evident. For instance, for a 
SMBH with M = 10^ AT© and a corotating torus with outer 
radius rout = 10^ AT and mass A/t = O.lAf (Aft = Af), 
the inclination with respect to the equatorial plane will de- 
crease, due to the hydrodynamic drag, for orbits with semi- 
latus rectum p > 35Af (p > 25Af ), while the EMRI signal 
will start being detectable by LISA already at a distance of 
w 50Af from the SMBH. Note, however, that unless one con- 
siders as the satellite an intermediate-mass black hole with 
m '^ lOOAf© around a 10^ Mq SMBH (a configuration which 
may be possible but about which too little is yet known), con- 
sidering EMRls at such large distances from the SMBH has 
the obvious drawback that the amplitude of the gravitational- 
wave signal will be proportionally smaller This will consid- 
erably reduce the detection volume, although the decrease in 
the event rate could be mitigated by the fact that weak-field 
EMRls are probably more numerous than strong-field EMRls, 
which are the ones accounted for in standard calculations of 
event rates. 

In general, we expect measurements of the evolution of the 
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inclination angle in the early phases of EMRIs to be a potential 
source of important information about the presence of thick 
tori which could not be detected by other techniques. More- 
over, because for any astrophysically plausible torus configu- 
rations the effect of the hydrodynamic drag becomes rapidly 
negligible in the very strong-field region of the SMBH (i.e., 
P ^ 5Af ), the presence of a torus would not prevent high- 
precision tests of the Kerr nature of the SMBH being per- 
formed. 

Although obtained with a simple model for the torus (i.e., 
with a constant specific angular momentum), the important 
feature that distinguishes the hydrodynamic drag from radi- 
ation reaction, namely the decrease of the inclination angle, 
cannot be affected by a change of the specific angular mo- 
mentum distribution (we recall that £ must be increasing with 
radius for stability). Such a feature, in fact, is simply due 
to the conservation of the momentum of the satellite in the 
6'-direction during accretion and to the dynamical friction of 
the fluid: both effects force the satellite to smaller inclina- 
tions by reducing its 6'-velocity. However, the calculation of 
the magnitude of the hydrodynamic drag and how it compares 
with radiation reaction for more general disk models is not 
straightforward. 

Tori built with increasing distributions of specific angular 
momentum, in fact, would have two substantial differences 
with respect to those considered here. Firstly, the separation 
between the specific angular momentum of the torus and the 
Keplerian specific angular momentum will generally decrease 
for orbits corotating with the torus, thus reducing the relative 
motion between the satellite and the fluid and consequently 
the hydrodynamic drag, whereas it will increase for orbits 
counter-rotating relative to the torus, thus enhancing the hy- 
drodynamic drag. The magnitude of this effect depends on the 
precise angular momentum distribution considered and rough 
estimates can be made assuming a power-law for the specific 
angular momentum, i.e., i/M ^ (r/AI)", with a < 0.5 for 
the torus to have an outer radius and a cusp IBTl [3911 . Using 
the general formulas reported in Sec. [Ill it is easy to check that 
for r between 20M and 50A/ and prograde orbits the relative 
motion decreases by ~ 20 — 30% for a = 0.1 and by at least 
95% for a — 0.4 (this significant decrease is due to the fact 
that for a = 0.4 the center moves to a radius reenter ^ 27 M, 
just in the middle of the radial interval which we are consid- 
ering). Conversely, in the same radial range the increase for 
counter-rotating orbits is of about 8 — 10% (30 — 45%), for 
a = 0.1 (0.4). Secondly, the density in the inner parts of 
the torus will generally decrease. Using again the expressions 
in Sec. |II] it is easy to check that the density decreases by 
about 18% (97%) at r - 20Af and by about 11% (90%) at 
r ~50M, fora = 0.1(0.4). 

Overall, therefore, the decrease of the inclination angle due 
to the hydrodynamic drag could be detectable by LISA also 
for non-constant I tori, especially if £ varies slowly with the 
radius or, if I varies rapidly with the radius, if the EMRI is 
counter-rotating relative to the torus. 

Finally, let us comment on two further effects that can in 
principle occur in the systems considered in this paper First, 
the motion of the satellite will be influenced by the gravita- 



tional attraction exterted by the torus. This is clearly a con- 
servative effect, and cannot therefore influence the infall of 
the satellite towards the SMBH, which is instead regulated by 
the dissipative forces (radiation reaction and hydrodynamic 
drag). However, this effect can in principle cause the perias- 
tron to advance, thus introducing a phase-shift in the gravita- 
tional waveforms. (Note that a similar advance is caused by 
the conservative part of the gravitational self-force [62].) To 
calculate the order of magnitude of this effect, let us consider 
for simplicity a thin disk of outer radius Tout, mass M-q and 
constant surface density, and a satellite of mass m located on 
the equatorial plane at a distance d <^ rout from the central 
SMBH, the mass of which we denote by M. The potential 
energy of the satellite due to the gravitational field of the disk 
can be easily calculated to be, up to a constant and to leading 
order. 



U : 






(70) 



This potential energy can be used to compute the Newtonian 
periastron precession of the satellite's orbit during a revolu- 
tion [use for instance eq. (1) of Ref. il63ll . chapter 3, exercise 
number 3]: 
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for almost circular orbits. Using this equation and the well- 
known Netwonian formula for the revolution period, it is easy 
to check that, for orbits relevant for LISA, the total phase-shift 
accumulated in 1 year is <C 27r as long as rout ^ 10** Af and 
M = 10^ - lO'^A^Q. Because LISA is not expected to detect 
phase lags less than 1 cycle over its lifetime'^, this periastron 
advance and the consequent phase-shift cannot be observed. 
On the other hand, for disks or tori with rout ~ lO'^Af , this 
effect could in principle be marginally visible by LISA (espe- 
cially if Af= lO^Af©). 

A second effect which could in principle affect EMRIs in 
the presence of a torus is the spin of the satellite black hole, 
which increases due to accretion of the torus material. The 
satellite's spin couples with the orbital angular momentum as 
well as with the spin of the SMBH, but its effect on the motion 
is negligible over a timescale of 1 year [64], unless it is close 
to its maximal value (in which case it might be marginally 
observable) 1, 64] . 
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